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ABSTRACT 

We reformulate the zero- dimensional hermitean one-matrix model as a (non- 
local) collective field theory, for finite N. The Jacobian arising by changing vari- 
ables from matrix eigenvalues to their density distribution is treated exactly. The 
semiclassical loop expansion turns out not to coincide with the (topological) ex- 
pansion, because the classical background has a non-trivial iV-dependence. We 
derive a simple integral equation for the classical eigenvalue density, which displays 
strong non-perturbative behavior around N = oo. This leads to IR singularities 
in the large- N expansion, but UV divergencies appear as well, despite remarkable 
cancellations among the Feynman diagrams. We evaluate the free energy at the 
two-loop level and discuss its regularization. A simple example serves to illustrate 
the problems and admits explicit comparison with orthogonal polynomial results. 
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1. Introduction 



Large- iV random matrix models have been with us since the 1950s [1], and 
enjoyed a renaissance 1989 due to newly-found applications to random surfaces, 2d 
gravity and string theory [2]. In the late 1970s matrix models encountered gauge 
theories, and its two major calculational frameworks were developed: semiclassical 
or collective field theory of matrix eigenvalues [3], and the technique of orthogonal 
polynomials [4] . The latter has driven the recent successes, by admitting the inves- 
tigation of the double-scaling limit and of the connection to integrable hierarchies. 
Given that the semiclassical approach is particularly simple in the planar limit, it 
is somewhat surprising that, except for the one-dimensional model, it has not been 
extended to the loop level. This paper is intended as an attempt in that direction. 

To be sure, there is a well-known obstacle to using the collective field in a ^ ex- 
pansion of the free energy, namely the exact evaluation of the Jacobian generated 
by the change of variables from matrix eigenvalues to their density distribution 
[5-8]. However, recent work of Karabali and Sakita [9] suggests a way of doing 
just that. In this paper we follow their approach and convert to the Fourier rep- 
resentation of the functional delta- function defining the eigenvalue density p{x). 
Its conjugate field, X(x), gets thereby introduced into an effective action Sn[p,X] 
which, though non-local and non-polynomial, takes on a fairly simple form. An 
equivalent result has been obtained by Ben-Menahem [10], through mediating the 
Coulomb repulsion of the eigenvalues by a 2d photon field. 

In section 2 we derive the effective action and its equation of motion, for 
any finite N. The full saddle-point equation contains an explicit term, which 
obstructs any attempt for an exact solution. At N = oo, it reduces to the airfoil 
equation, which admits exact solutions vanishing outside a compact support a 
[11,12,3]. For N<oo, however, strong non-perturbative (in -^) behavior appears 
off a. Section 3 presents the semi-classical loop expansion around the classical 
background at finite N, i.e. one-loop determinant, vertices and propagators. The 
Feynman graphs for the free energy Fn are drawn at the two-loop level. The exact 
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but formal expressions are then expanded around iV = oo, in section 4. At first 
we hide the implicit ^-dependence of the saddle-point density p(x) and obtain 
remarkably simple results for the loop corrections to F^, due to large cancellations 
among the Feynman diagrams. We then develop p(x) in jj as well by iterating a 
system of airfoil equations. The expansion is beset by IR and UV singularities 
whose regularization is briefly discussed. Section 5 works out in detail a simple 
example and compares to results obtained by orthogonal polynomial techniques. 
The conclusions form section 6, and some technical material as well as a useful 
compendium of formulae is collected in three appendices. 



2. The Classical Eigenvalue Density 

We begin by formulating a collective field theory for the zero-dimensional her- 
mitean one-matrix model at finite N. Our starting point is the partition function 



Z N oc jd N Me- NtlV ^ 



(2.1) 



for an N x N random hermitean matrix ensemble, in a potential V. Upon diago- 
nalization M = diag(xi) this reduces to 



Zn '■= 



N r i r 

J] / dxi exp j -N J2 V(xi) + ln( 

i=l J ^ i i<j 



X{ Xj ) 



:2.2) 



We like to change variables from the matrix eigenvalues X{ to their density distri- 
bution* 



N 



P( x ) ■= jr^2 S ( x ~ x i) ■ 



(2.3) 



★ Obviously, this transformation is N : Hi but then p vanishes almost everywhere. 
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More precisely, we insert 



1 = fvp Y[6(p(x) - ^(x-xi)) 

= J J VpVX exp{* J dxX(x) [p(x) - ^Ax-Xi)} } 



(2.4) 



into (2.2) and express the action in terms of the density, 

S° N \P] = N 2 jdxp{x)V(x)-\N 2 jjdxdy p{x)f{x-y)p{y) + \N f{$) . (2.5) 

The self-interaction had to be regulated by replacing \nz 2 — > f(z) in eq. (2.2), 
choosing some suitable, i.e. symmetric and bounded, function /. 

Following ref. 9 we are able to perform the integration over Xi, 



Zn 



VpVX e - s * [p]+l f Xp 



N 



■i=l 



1 1 / dx 



dx e N 



7v A ( a; ) 



i N 



(2.6) 



and arrive at an effective action 



S N [p,X] = S° N [p]-i J dxX{x)p{x)-N\n J dx 



(2.7) 



which is not only nonlocal in the two real fields p and A but also non-polynomial in 
the latter. Interestingly, the constant mode of A can be integrated out exactly to 
yield the constraint 5(J p — 1) that was apparent already from the definition (2.3). 
However, we shall keep those modes in the measure for the time being. In principle, 
another constraint arises from the positivity of p. Perturbation theory about a 
strictly positive p, however, is insensitive to this restriction, and we will therefore 
ignore it in the fluctuations. 
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Our goal is to carry out a systematic semiclassical analysis of this peculiar one- 
dimensional field theory. To leading order in h we must determine the saddle-point 
configurations (p, A), where the action (2.7) is stationary. The first variations yield 



= N 2 V(x)-N 2 dyf(x-y)p(y)-i\(x) 

(2-8) 

= p(x) - e~¥ X{ - x) /fdye-V X ^ 

where p comes out to be properly normalized, J p = 1. The second equation 
determines A up to a constant, 

X(x) = iN In p(x) + A , (2.9) 

a result which may be inserted into the first equation* Subtracting the value at 
x = to get rid of the constant and deregulating the Coulomb repulsion we get^ 

[V(x)-V(0j\- Jdy[\n(x^) 2 -\ny 2 ]p(y) + ±[\np(x)-\np(0)] = . (2.10) 

Together with normalization and positivity, this equation describes the classical 
eigenvalue density for any finite N. It is noteworthy that (2.10) is not homogeneous 
in jy-, so its solution cannot be, either. The saddle-point equation is more easily 
studied in its differentiated version, 

\V{x) + ^6*) = /-^-p(y) = nH x [ P ] (2.11) 
p j x y 

where -f denotes Cauchy's principal value of the integral. The r.h.s. is known as 
the Hilbert transform TL (of p) which has been thoroughly investigated [13] (see 

* The first equation subsequently fixes the constant Ao for a given solution p. 
f See also rcf. [10]. 
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appendix A for details). At this point let us merely state that 7i is uniquely 
invertible on the class L p , p > 1, by 



Kx 1 



■--f-^-m ■ (2.i2) 

7r J x-y 



At iV = oo our equation (2.11) has been widely studied and solved [11,12,3], 
and it was learned [14] that a unique solution extends to iV < oo. Unfortunately, 
the equation is not easily solved for finite N. Even in the large- N limit some care 
is required, as shown by the following. For instance, should we drop the p'/p-term 
since it is down by -^? A little inspection reveals that such a step is in general not 
consistent with the asymptotic large- \x\ behavior of the equation, which demands 
a i fall-off for the l.h.s. due to the normalization of p. In fact, the formal iterations 

Kx) = \\V{y) + ^{y)] 

tx 1 J x-y l z p J ^ 2 13 ^ 

p(x) ~ exp{-N[V(x) - f\n(x-y) 2 p(y)]} N ^°° x 2N e~ NV ^ 

show all terms in (2.11) to be of the same order for \x\ ^> 1 (unless V ~ hire 2 ). The 
situation is different, however, near the minimum of the potential where most of p 
is concentrated? When p is 0(1) the p'/p-term may safely be neglected or treated 
as a ^ perturbation in (2.11). To the same approximation one should also drop 
the exponential tail of p under the integral. At iV = oo this amounts to setting p 
identically to zero beyond a certain distance from the minimum of V. The cross- 
over points x = a,b between the exterior region (where p is 0(e~ N )) and the 
interior domain a = (a, b) (where it is 0(1)) are determined by the normalization 
of p; they occur where H[p] can no longer match the growth of V' . For iV — > oo, 
p — > po, it is therefore indeed consistent to simplify eq. (2.11) to [3] 

b 

-^—p (y) = 7iT x [p ] forxe(a,b) 

a 

po(x) = for x ^ (a, b) 



| For simplicity we restrict ourselves to single-well potentials. See, however, ref. 14. 
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This relation is known as the airfoil equation and determines the vorticity po related 
to a given velocity field V' along the airfoil a (b—a is the span of the wings) [11]. 
It is solved by inverting the finite Hilbert transform % 



p (x) 



y/(b -x)(x- a) 
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dy It// 



x-y 



V'{y) y/(b-y)(y-a) 



(2.15) 



where the arbitrary constant 



7 = \ J d v po(v) 



(2.16) 



arises due to the one-dimensional kernel 



\, ; of ^"and is set to 7 = ± by 

y/(b-x)(x-a) n 

the normalization of po- The inversion formula (2.15) is guaranteed to produce a 
solution to (2.14) as long as V' G L^ +e ; for V' G L 1+e we also need po G L 1+e 
in (2.15) [11, 12]. In the case of polynomial potentials the resulting density po(x) 
is seen to vanish at the endpoints of a generically as (x — a) 1 / 2 and (b — re) 1 / 2 . For 
more on the finite Hilbert transform we refer to appendix A. 

Having solved the saddle-point equation in the large- N limit, one is tempted 
to try and expand the complete solution p of (2.11) around po- A -^--expansion 
appears logical on a but is bound to fail in the exterior region since e~ N has 
an essential singularity at = 00. Instead, a series in powers of e~ NV seems 
more appropriate there. Either type of expansion, however, must break down at 
the cross-over points, because a smooth function p is to develop from one that is 
singular there. For finite N the transition from 0(1) to 0(e~ N ) occurs in a finite 
intermediate region of size A~ 2 / 3 by a generalized Airy function of the scaling 
variable {x— 6) A 2 / 3 . Still, one might hope that a -^--expansion of the classical 
action needs p only on a, due to its non-perturbative nature off a. We shall pursue 
this line of argument in Sect. 4. 
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It remains to evaluate the classical action in terms of p. Employing eq. (2.9) 

_i_ c 

as well as In J e N = — -^Aq we find, in agreement with ref. [10], that 



S N \p,X] = N 2 JpV-^N 2 JJpfp + N Jp\np+\Nf{0) (2.17) 

which may be rewritten as 

N 2 f N 2 N f N N 

Sn[p,\] = -y Jp(V-f) + ^V(0) + - yplnp + -lnp(0) + -/(0) 



= ^ IJ p(x)[f(x^y)-f(x)-f(y)}p(y) + N\np(0) + ^m 

(2.18) 

by making use of the saddle-point equation (2.10). Given the non-trivial N- 
dependence of p it becomes clear that the semiclassical expansion will not cor- 
respond exactly to an expansion in On the other hand, it is encouraging to 
note that potentially divergent pieces like J In p have dropped out. The eventually 
singular /(0) term will get cancelled in the following section. 

3. Quantum Fluctuations 

We now set out to develop the loop expansion of the theory defined by eqs. (2.5) 
and (2.7) in the standard fashion. We eliminate the constant modes from the 
measure and split 

p(x) = p(x)+p\x) X(x) = \{x) + \'{x) (3.1) 

with fp' = = f\'. Since Sn is quadratic in p its fluctuations are gaussian and 
affect only the propagators. Hence, only vertices for A' are present. Let us write 

Z N = Z d Z!_ loop Z+ Z cl = e ~ SN ^ Z!_ loop = [det%]"5 

(3.2) 

where In Z + equals the connected Feynman diagrams, and det' reminds us to omit 
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the constant modes. The fluctuation operator S'^ at the saddle point reads 

f-N 2 f(x-y) -iS(x-y) \ 

Sn(.x, V ) = 1 r . (3.3) 

y -t5(x-y) [p(x) 5(x-y) - p{x)p{y)\ J 

and acts on {p'(y), X'(y)}. The simple form of the off-diagonal blocks of S'^- facili- 
tates the calculation of the determinant* 



Zi-ioop = det'[5(x-y) - Nf(x-y)p(y) + Njdz f(x-z)p(z)p(y) 

(™N k r nfc 
= exp<^ 2^ — tT\f(x-u)p(y) - Jdz f(x-z)p(z)p(y) 

^k=i 



(3.4) 



where we have formally expanded in powers of N (not -^). In the obvious graphical 
representation this gives rise not only to the usual loop diagrams but also to loops 
"cut open" at any line, due to the second term in the trace. Each term in the 
above expansion seems to be well-behaved except for the very first one, ^tr/p = 
y/(0)j which in fact just cancels the eventually singular term in the classical 
action, eq. (2.17). Of course, one would rather like to expand in powers of but 
let us defer this to the next section. 

Instead, we complete the construction of the Feynman rules. To obtain the A' 
propagator we need to invert S'^, 



X'(x)X'(y) = N[(p(x)5(x-y)-p(x)p(y))-^f- l (x,y)] 



-l 



= -N 2 f(x, y) - N z jj f(x, z) { P {z)5{z-w) - p(z) p(w)) f (w , y) - . . . 

(3.5) 

again formally expanding in N. The corresponding graphs are the same as in (3.4) 
before taking the trace. 



* det( A a g)~ dot [1 + A ■ B}. Since 1 is a zero mode of B we may omit the primes, 
f The inversion of / will be discussed in the next section. 
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Turning to the vertices, the A' expansion simplifies at the saddle point as 



AS = -Nln Jdxe-^ X{x) = iX - N In Jdx p(x) e"I A » 



r=l 



■ LXJ 1 



(3.6) 



with the help of e w / j e iv = p. The constant mode of A' appears only at 
0(1), as it should. Note that this is a genuine jj expansion because in powers 
of -^A'. The (9(1) and O(j^) terms have already been taken care of in (2.8) 
and (3.3). 

It is instructive to compute the two-loop correction to the free energy 



F N = -\nZ N = S N [p, X] + ±tr'\nS%- In Z+ . (3.7) 



Writing AS = N N k Sk, with k lines emanating from Sk, we have 



-S 3 /N 2 -S 4 /N 3 -. 



exp{-iV- 3 (S 4 ) + \n- a (S 3 S 3 ) - . . .} (3.8) 



where (. . .) denotes the free field average performed with the propagator (3.5). 
With a little combinatorics Sk can be read off eq. (3.6), e.g. 



3 r 



^3 = 

s 4 = 



H) 

3! 

(-0 4 



-fp\* + 3jpX n Jp\' - 2(J P X 



12 1 



i/\3 



4! 



-JpX' 4 + 4fpX«fpX' + 3(/pA' 2 ) 2 - 12 fpA' 2 (/pA') 2 + 6(/pA') 4 

(3.9) 

which are depicted in figure 1. Obviously, the parts of Sk are labelled by the 
partitions of the number k into integers. The A' contractions produce 21 diagrams 
which are displayed in figure 2. Due to the non-locality of the vertices, the terms 
"loop" and "vertex" take on two distinct meanings, namely one counting powers 
of h and denoting Sk, and the other describing individual diagrams. In order to 
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avoid notational confusion, we distinguish between the two and refer to figure 2 
as the two-loop contributions built from vertices Sk, represented by 6 two-loop, 9 
one-loop, and 6 tree-level diagrams consisting of 2 or 3 lines and up to 6 nodes. 
Note also that all contributions to \nZ + are connected whereas the constituting 
graphs may not be. It appears that we have derived a sensible set of Feynman 
rules (3.4)-(3.6) for our non-local theory. 

The functional AS in (3.6) has the property that an ^-independent shift of 
A' leaves the vertices Sk form-invariant. We may exploit this fact to simplify the 
structure of Sk, by changing the field variable A' — > A := A' — f p\'* Since JpX — 0, 
a large number of terms disappear when Sk are expressed in terms of A. We are 
left with the parts corresponding to partitions of k into integers > 2 only. The 
21 diagrams of figure 2 compactify to the five graphs without univalent nodes. 
In particular, tree graphs can no longer be present. The price we pay for this 
simplification is a more involved propagator for A, which actually consists of four 
terms when rewritten in terms of the A' propagator. 

We gain, however, an argument showing the cancellation of all tadpole infini- 
ties. The only potentially divergent diagrams arise from self-contractions, 

X'(x)X'(x) = -N 2 f(0) + finite terms . (3.10) 



Any A-propagator ending in a univalent node with a weight of constant x p is 
seen to vanish, which eliminates all /(0) terms from tadpoles. The singularities 
arising from self-contracted multiply valent nodes appear to cancel systematically 
between different graphs (see figure 2 for example). This argument holds true for 
any constant part of X'(x) A' (a;), whether we expand in N or not. 

Expanding the full propagators in the loop graphs a la (3.5) leads to an in- 
finite series in powers of N, which ruins the expansion from (3.6). Indeed, 

* The Jacobian of this transformation is compensated by the change in the one-loop determinant. 
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absorbing the leading iV-dependence of X'(x) A/ (y) into the vertices in (3.8) re- 
places Sk/N^ 1 — > NSk, for any k, so that all terms in AS* contribute at lowest 
order in N. It therefore appears that we need to sum the loop expansion in order to 
consistently develop Fjv in powers of N. On the other hand, Fjy has a logarithmic 
singularity at iV = (at least for a quadratic matrix potential) which makes an N 
expansion ill-defined anyway. 

Nevertheless, the semiclassical expansion seems to be well-behaved for finite N. 
It is merely dangerous to develop the individual loop contributions in N (or as 
we shall see soon). Unfortunately, for practical purposes we need to compute the 
determinant (3.4) and the propagator (3.5) somewhere, and this seems feasible 
only at the extrema N = or N = oo. Another option is to put the system in a 
box (in eigenvalue space) of length L and estimate determinant and propagator in 
Fourier space. We shall come back to this approach later on. 

4. ^ Expansion 

This section is devoted to an attempt to set up a systematic expansion of 
the free energy (3.7), using the previous results (2.11)-(2.18) and (3.4)-(3.6). We 
will proceed in two stages. First, we arrange according to explicit powers of and 
hide the implicit iV-dependence of p. Second, we expand p as well and insert the 
result into the coefficients obtained earlier. 

The classical part of the free energy given by eqs. (2.17) or (2.18), it remains to 
develop the fluctuation determinant (3.4) and the full A' propagator (3.5) around 
iV = oo. To this end one needs to invert NS'^\\\ = p(x)5(x—y) — p(x)p(y) away 
from the constant mode. As expected, this operator has a constant zero mode 
since Jp — 1. Hence, its inverse orthogonal to that mode is simply 

[p{x)5{x-y) - P (x)p(y)]^ = ^ 5(x - y) . (4.1) 

f It is convenient to make a change of basis, A' = £,/\fp, so that the kernel of S(x—y) — 
\J p{x)p{y) is spanned by \/~p. Because y ' p{x)p(y) projects onto the kernel, our operator 
projects onto the complement and is inverted by the identity there. 
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On the level of the determinant this can also be seen in Fourier space, from the 
identity 

det '[pk+i - pk Pi] = det[p k+ i] for p = 1 . (4.2) 
The propagator (3.5) is then easily written out as 

X'(x)X'(y) = N^5{x-y) + j^f^y)^ + i/)/- 1 )/" 1 ) + O(^) ■ 



(4.3) 

The one-loop determinant (3.4) works out analogously 



_ i 

Zi-loop — det 



-Nf(x-y)p(y)] 5 x det [s(x-y) - ^f-\ X: y) 

oc expj-^lniV- ^ ln p| x ex p|Xl tr P^)" 1 /"^ 2/) j 

(4.4) 

where rj is a regulating constant, and we dropped a singular constant determinant. 
Together with eqs. (2.18) and (3.6) this yields a systematic but formal expansion of 
the free energy (3.7). Each A;- loop contribution starts at 0(N~ k+1 ) and is made of 
k- and lesser-loop graphs, many of which are disconnected. The ensuing Feynman 
rules are given in figure 3, but the combinatorical coefficients and the occurrence 
of the nodes have to be read off the expansions (3.6), (4.3) and (4.4). Clearly, 
the expansion is not just in but also in which will be responsible for IR 
divergencies later on. 

We have yet to give a less formal meaning to the "propagator" Since the 
convolution with -^d x f(x,y) — > ^P^z^ is just the Hilbert transform (2.11), the 
deregulated f(x—y) — > \n(x— y) 2 is easily inverted from (2.12)J 

r\x,y) - -_L P ^_ * = "^P— (4-5) 
27r z x — y ay 2ir A ax x — y 

on the space 



$ The symbol P denotes the principal value, J P . . . = . ., or P^ = j — in5(x). 
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Employing the results of appendix A we are now able to compute the one- 
loop traces of eq. (4.4). The derivatives get distributed over the graph, and, after 
evaluating 



Jj^^^-^^y) = /Wi] & - *(o) y 0i 02 



P ^(x)— 2 (y)— 3 (Z) = /(H / [0l]02-01^ / [02])H[0 3 ] 

x-y y-z J 



z — X 



(4.6) 



and so on, we discover that tr [p -1 / -1 ] fe vanish upon symmetrization over cyclic 
permutations! ' } Hence, the one-loop correction to Fn comes from the first (singu- 
lar) factor of (4.4) only and is given by 

itr'ln^ = \r]\nN + \ Jdxhip(x) . (4.7) 

By the same calculation, self-contractions \'(x) V (rr) do not receive corrections 
beyond the leading term. 

To order iV _1 we have to review our two-loop calculation of the previous sec- 
tion. It suffices to replace all full propagators of figure 2 by N p(x)~ 1 5(x—y). In 
terms of the Feynman rules of figure 3 this means that we compose our diagrams 
with dashed lines only. Thus any propagator loop yields a singular factor of 5(0), 
leaving us with connected graph components proportional to 5(0) 2 p _1 , 5(0), or p. 
It is quite remarkable that all but one of the divergent graphs cancel between (S4) 
and (SsSs)^ However, the single genuine two-loop diagram survives (together 
with some finite contributions), 

-\nZ. 



J + 



12 y ' J p{x) 12 



A7 + • ^ 



It is then quite clear that a fc-loop contribution to In Z + will yield (products of) 
expressions of the form 5(0) p J p 1_p , with < p < k. However, all Feynman graphs 



§ We have checked this for k < 3. 

The propagator self-contraction cancellation argument of the previous section does not 
even apply since ^'(x) ^ (x) = Nd(Q)/p(x) is not constant. 
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appear to cancel each other except for those with p = and p — k. The latter are 
simply given by all genuine (connected) /c-loop graphs, presumably without self- 
contractions. 

It should not have escaped the attentive reader that the expressions (4.7) 
and (4.8) are divergent. In fact, we encounter two different types of singulari- 
ties. First, the regulator rj as well as the factors of 5(0) arose from evaluating trl 
and thus count the (infinite) number of degrees of freedom. Second, the integrals in 
both expressions are singular since p ~ e~ NV for \x\ — > oo induces J hip ~ iV J V 
and J p 1 ^ J e +NV . The latter infinities require an IR cutoff; most convenient 
is a restriction of the eigenvalue range to a finite interval, say x e [ — "with 
open boundary conditions. Putting the system into a box does not, however, cure 
the UV singularities mentioned first; those ask for a mode cutoff, say \k\ < M, for 
some large M. A natural choice for the value of M would seem M = aN+(3. We 
have estimated the resulting leading behavior of the one-loop determinant (3.4), 
with some details given in appendix C. The upshot is 



itr'ln^ 



reg 



[(a+l)ln(a+l) - a\na]N + Q(\nN) . (4.9) 



The O(N) term originates from ^lndet[— /] which we dropped in (4.4). Since 
we did not probe the p dependence, L does not appear in our result. Further 
numerical analysis is required to elucidate the cut-off dependence of the regulated 
loop corrections. 

The coefficients of the series we have obtained are (singular) functionals 
of the saddle-point density p. To complete the jj expansion we finally develop p 
about N = oo as well, thus blinding ourselves to its "non-perturbative" parts. In 
perturbation we set p(x) = for x ^ a and expand 

oo 

p(x) = ^ iV~* pi(x) ioixea . (4.10) 

i=0 



* We have checked this also on the 3-loop level, with an outcome of j^5(0) 3 J{pN) 2 + finite. 
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Next, we define a sequence of potentials 



v a = v , y 1 = in P „ , £jv-'y j+1 = £fcL_ \Y.N-e£\ 

i=l r=l ^s=l ' 

(4.11) 

by expanding In p. The saddle-point equation (2.11) then turns into the iteration 



b 



i V f(x) = -f-^- Pi(y) for x G (a, b) (4.12) 
J x-y 

a 

since V% is a function of {po,Pb • • • , Pi-i}- The constant part of pi is fixed by the 
normalization conditions jpi = <5i,o- Clearly, the iteration consists of repeated use 
of eqn. (2.15) for T _1 , with 7 = after the first step. 

Having computed the coefficients in (4.10) to a desired order in we are able 
to unravel the complete large- N expansion of the free energy. First, the above 
equations can be employed to simplify the jj series of the classical action (2.18): 

Sn[p,\} = ^Jp Q [V(x)-V(0)-f(x)]+^f(0) + Njp \np + ^ jp 1 \np + ... 

= \N 2 ff p ( X )[f( X -y)-f( X )-f(y)]p (y) + lNf(0)+N\np O (0) 
+ NJJ p (x) [f( X -y)-f(y)] pi (y) + \ JJ pi (x)f(x-y) Pl (y) + 0(jj) . 

(4.13) 

Second, we insert the expansion (4.10) into the expressions (4.7) and (4.8) for the 
loop corrections. However, by setting p ~ e~ NV to zero off a we are dropping 
infinite contributions ~ NjV and ~ Je +NV to the integrals of (4.7) and (4.8), 
respectively, of which the first is perfectly perturbative. 

At this stage it is quite obvious that a brute-force expansion around iV = 
00 makes no sense. In view of the interference with strong e~ N behavior one 
must keep at least the iV-dependence of p intact. Only then may one hope to 
establish a large- N series for the free energy, provided the regulation problems can 
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be overcome. Under these circumstances it is hard to see how a serious comparison 
with orthogonal polynomial results could be achieved beyond the leading order. 



5. An Example 

We have tried to develop the semiclassical expansion for the zero-dimensional 
hermitean matrix model with an arbitrary potential V. In order to explicitly 
demonstrate the difficulties of this scheme we now become more specific and choose 
a simple potential, 

V(x) = \x 2 + gx A , g>0 , (5.1) 

for which we can evaluate the free energy Fjv by the method of orthogonal poly- 
nomials. From the literature [4] we extract the value for F^(g)—. F/v(0), and add to 
it the result for the quadratic potential (g = 0), as calculated in appendix B. The 
first orders in yield 



F N = 



jj + ^(a 2 -l)(9-a 2 )-^hia 2 



iV 2 -N In N 

(5.2) 

+ [l — In 2ix\ iV — — In iV + c + — ln(2 - a 2 ) + 0{N- 2 ) 



with 

12#a 4 + a 2 -l=0 « a 2 = (-1 + y/l + A8g)/2Ag (5.3) 

and Co ~ —0.7535, a universal constant. 

Let us see how much of this the collective field is able to reproduce. The 
iterative solution of eqs. (4.11) and (4.12) proceeds along standard lines. To lowest 
order eq. (2.15) produces [3] 

po ( x ) = ^ (1 + 8ga 2 + 4gx 2 ) V '4a 2 - x 2 0(4a 2 - x 2 ) , (5.4) 

a polynomial deformation of Wigner's semicircle law. The next step in the itera- 



te 



tion (4.11) and (4.12) yields 

2x 



V((x) = 



x 



x 2, + n z 4a z — x l 



with n 2 = lif^ = f±^a 2 (5.5) 



which is to be plugged into (2.15) to find the correction pi = p^+p^\ With the 
help of appendix A* the first term gives 



1 



7T V4a 2 — x" 1 



A/(4a 2 + n 2 ) n 2 



6(4a 2 - x 2 ) 



(5.6) 



x 2 + n 2 

but the second term is not integrable at x — ±2a; stretching the meaning of (4.12) 
we might set 

1 0(4a 2 -:r 2 ) 



pf\ x ) = - [5(x - 2a) + 5{x + 2a)] - 



(5.7) 



4 L ~ V ~~ _/ ~ v " " n br V4a 2 - x 2 

These results continue to hold true for — ^ < g < 0; in this regime 2 > a 2 > 1 
and —8 > n 2 > — oo. At the lower bound \n\ — > 2a, and the well-known additional 
zero develops at the edge of a. Evidently, the iteration (4.12) leads to ever worse 
singularities at x — ±2a. This merely reflects the breakdown of the expansion 
at the edges and shows that our neglect of the non-perturbative tail of p was 
illegitimate. 

Nevertheless, we insert the above findings into the functionals (4.13), (4.7) 
and (4.8). The classical action becomes 

S N [p,\] = N 2 [^ + ±(a 2 -l)(9-a 2 )-Una 2 ] + \Nf(0) 



+ N 



- In 7i + a 2 (i-ln 2) + In + (l-a 2 )(^-ln 2) (5.8) 
+2 In 1(1+6) + i(2V){^ + \{l^ 2 ){Uf] + 



with b 2 = 3|rp^ and a singular outcome for 0(1). The one- loop correction yields 



itr'ln^ = ir/lnA^-2a[l + lnf -ln 2 ^ -2( C - 1 arctanc- 1) 



(5.9) 



where c 2 = 2 \^ a i , which does not carry much meaning in view of the singular 



* It it useful to remember that - i^] dx 

7T J— 1 



a? 2 ±a 2 



-l + e^-^viia- 2 . 
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classical contribution at the same order. 

The comparison of these semiclassical "results" with the exact formula (5.2) 
is poor. Although the leading terms agree (as is well-known) we have failed to 
reproduce the subleading contributions. The logarithmic iV-dependence may still 
be restored by taking ln(x-^/) 2 = In N~ 2 when \x-y\ < iV -1 , and by setting r] = — |. 
The O(N) term is (incorrectly) a-dependent and even at g — (a 2 = 1) does not 
quite match up. We do, however, expect corrections of 0(N) and below from a 
proper regularization of the one-loop determinant (see (4.9)), hence only the first 
line of eq. (5.8) can be trusted anyway. Presently, the first non-trivial subleading 
term of F/v(g)— F/v(0) cannot be probed due to the divergencies, although it is 
interesting to note that the logarithmic singularity at a 2 = 2 already shows up 
correctly. 



6. Conclusions 



We have formulated the collective field theory for hermitean matrix models 
away from the large- N limit. By treating the Jacobian from the change of variables 
from eigenvalues X{ to their density p(x) exactly, we arrived at a particular non- 
local and non-polynomial theory of two real scalar fields, p and A. Nevertheless, 
we were able to set up a loop expansion for the free energy, at any N < oo. The 
classical eigenvalue density, p, is determined by a simple integral equation whose 
large- N limit becomes the usual airfoil equation. Its action naturally contains the 
entropy term for the Dyson gas of eigenvalues. The H expansion about the classical 
background yields sensible Feynman rules, including somewhat formal expressions 
for the one-loop determinant and the full propagator. It has to be stressed that the 
semiclassical (loop) expansion is not identical to the topological (-^) expansion, 
because the "equation of motion" for p is not homogeneous in N, thanks to the 
Jacobian. We studied the collective field theory in some detail up to the three-loop 
level. 
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The attempt to develop our non-local field theory around N = oo met with 
three difficulties. First, the classical density p(x) can be expanded in only near 
the minimum of the matrix potential V(x); elsewhere it behaves asymptotically 
as e ~ NV ( x \ a strong non-perturbative (in jj) signal. The mixing of both regions 
due to the defining integral equation destroys the expansion everywhere; only the 
leading term p — > po is non-singular. Second, UV singularities oc (trl) fc appear at 
k loops since the propagator is a delta-function to leading order. Amazingly, most 
of the divergent graphs cancel each other, except for the connected diagrams with 
the maximal number of loops and no self-contractions. As a result, the free energy 
is not just expanded in but also in Third, the ensuing integrals J p 1 ~ k do 
not converge due to the exponentially small tail of p, producing IR divergencies. 
A simple recipe to deal with these problems is to restrict the eigenvalue space to a 
finite interval and to implement a high-frequency cutoff. Preliminary answers were 
obtained here at the one-loop level. 

Finally, a simple example was chosen to substantiate our abstract discussion 
and to compare with known results obtained by the method of orthogonal polyno- 
mials. As anticipated, only the planar contribution was reproduced entirely, but 
some features of the torus correction appeared as well. 

Although we focused on the zero-dimensional hermitean one-matrix ensemble, 
the discussions in this paper should carry over mutatis mutandis to multi-matrix 
models, the one-dimensional theory, and to unitary matrix models. 

It must be said that we did not yet achieve our goal of deriving a practical 
algorithm which systematically computes corrections in the semiclassical scheme. 
However, the necessary steps have been taken, and the details of our failure are 
no more puzzling as they are intriguing. We certainly hope to gain further insight 
into this outstanding problem in the future. 

Acknowledgements: I wish to thank D. Karabali, J. Schiff and S. Wadia for useful 
discussions. 
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APPENDIX A 



We give the basic notions of the infinite and finite Hilbert transforms, including 
some useful properties. Hilbert was the first to notice the (skew-)reciprocity of 

oo oo 

u(x) = -J-^—v(y) , v(x) = ---f-^—u(y) ; (A.l) 
vr J x-y 7T J x-y 

— oo — oo 

u is said to be conjugate to v. The (infinite) Hilbert transform 

oo 

u(x) = H x [v] := - J -^-v(y) (A.2) 
7i J x-y 

— oo 

is related to Fourier integrals and to the Laplace transform and received an ex- 
tensive treatment in Titchmarsh's book [13]. For v £ L 1+e , Tt[v] exists almost 
everywhere and also belongs to L 1+e , so its inverse is 

v(x) = H- l [u\ = -H x [u\ . (A.3) 

If v G L, H[v] still exists almost everywhere but does not necessarily belong to L. 
We may, however, set = 0. Furthermore, it is easily seen that 

H x [v(ayH>)] = ^ H ax+b [v(y)] . (A.4) 
Some interesting properties of Ti generalize to the finite Hilbert transform, 

it J x-y 

a 

which plays an important role in aerodynamics and has been investigated by a 
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number of authors [11, 12]. Most useful of those is Tricomi's convolution theorem, 



T{viT[v 2 ] +v 2 T[ Vl ]} = 7>i]7> 2 ] -viV2 , (A.6) 



valid for V{ G L Pi with Pi > 1 and ^ + ^ < 1. It is equivalent to the interchange 
formula 



b b b b 



a 1a v i(v) v 2( z ) _[, J, vi(y)v 2 {z) 2 

dy4dz- = -fdz-jdy- - ir v 1 {x)v 2 {x) (A.7) 



of Hardy and Poincare (more generally, v\(y) v 2 (z) — >• F(x, y, z)). Under the same 
conditions we also have Parseval's formula 

[dxJdy Vl{x)V2{y) = b [dyJdx Vl{x)V2{y) (A.8) 
J J x-y J J x-y 



which may also be written as 

b 

Jdx (v 1 r[v 2 ]+v 2 r[vi]j =0 , (A.9) 

a 

demonstrating the skew-adjoinedness of the T(and H) operations with respect to 
the standard scalar product. Another convenient relation is 

r x [v] = T x [v'\ (A.10) 

which allows us to commute a derivative past the transform, as in (4.5). All these 
relations, (A.6)-(A.10), hold for H as well. 



21 



Conversely, the linearity property (A. 4) cannot hold for T; rather 

x — b 



Till] 



7T 



rr — a 



(A.11) 



so that (A. 6) with v\ — l yields T{T[t>]} ^ —v, in contradistinction to Ti. This 
slightly complicates the inversion of T. However, it is easily calculated that 



%[ 



y/(b-y)(v-a) 







T x [y/(b-y)(y-a)\ = X - 



a+b 
2 



(A.12) 



on the interval x G (a, 6). With these data a clever application of the convolution 
theorem (A. 6) then establishes [12] 



y/(b -x)(x- a) 



7 



hi 



dy 
x-y 



u 



(y) V( b -y)(y~ a ) 



(A.13) 



where, using (A. 8) and (A.12), 



7 





= ~ J dx T x\ u \ 



(A.14) 



is an arbitrary normalization constant. The conditions of the convolution theorem 

4 

restrict the validity of the inversion (A.13) to u e L3 +e , but the uniqueness of the 
kernel (A.12) holds in L 1+e . For u E L 1+e the r.h.s. may not be in the same class; 
on the other side, any inverse of Tin L 1+e must be of the form (A.13). Furthermore, 
given u G u] is necessarily orthogonal to the zero mode, 





J dx 



y/(b-x)(x-a) 



%[u] 



(A.15) 



by virtue of the skew-adjoinedness (A. 9). This allows for alternative representa- 



-p i 

* A note of caution: outside (a, b) the r.h.s. acquire an additional term of ±e[(x-b)(x-a)] 2 , 
respectively, with e = 0(x— b) — Q(a—x) = ±1. 
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tions, e.g. 



b 



T-\u] = - 7 - i /-*L «(„) 4^)^) (A.16) 

^(fe-xJCx-a) Ws-y y/(b-y)(y-a) 

which also goes to (A. 3) when (a, b) — > (—00, 00). 

Unfortunately, our nonlinear extension of the Hilbert transform, 

u(x) = H x [v]-P^(x) (A.17) 
does not seem to have been investigated in the literature. 

APPENDIX B 

In the case of a gaussian ensemble (V = \x 2 ) of hermitean random matrices 
the free energy Fjv for any A can be found exactly, with orthogonal polynomial 
techniques. As usual [4] 



e 



N r 1 - _ N_1 

- f » - n e-^s^? n^-^-) 2 = ^ (b.i 



i<j fc=0 



where h k is the norm of the fcth monic orthogonal polynomial, P k {x) = x k + . . .. For 
our quadratic potential the P k are just proportional to the Hermite polynomials, 

P k (x) = V2N~ k H k (x^N/2) , (B.2) 

whose norms are well-known, 

h k = f+™dx e-* Nx * P k (x) 2 = N~ k -^k\V2^ . (B.3) 

Hence, we find 

N 

e -F N = N -N*/2 {2n) N/2 "Q fc , (B 4) 

k=l 

which is the exact result. 
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We go on to expand 



F N = ±N 2 \nN - t^N In 2ir - J2i (N-k+1) In k (B.5) 
in powers of -fa, with the help of Stirling's asymptotic formula 



lnr(z) = (z-I)ln^-^ + Iln2vr + Ei 1 ^1) * 1-2 * + Ra(z) 



(B.6) 



2s(2s-l) \z\ 2s ~ l cos^-^iargz) 

for large values of \z\ as well as 



N N+1 

J^klnk = fdt]nT(t) + l(N+l) (N-]n2n) (B.7) 
k=i { 

p+1 

Jdt\nT{t) = plnp - p + iln27T . (B.8) 



p 

T2- 



The N In N terms cancel, and after collecting all the pieces we arrive at the 
asymptotic series 



F N * -N> - N In N + (1-ln 2«)N - - In N + q, - g 3^)^-2) < M > 
where 

Co ~ — | ln27r — ^ + J]^ 1 

+ |P 2 - |(p 2 ^) Inp + ET 2.(5^^1-2) (B.10) 
= -0.7535173895042... 

is a universal constant, defined ever more accurately for ever larger values of the 
cut-off parameter p. As expected on general grounds, this is an expansion in 
except for a few terms at low orders. 



* B 2 k denotes the Bernoulli numbers: B 2 = ^, B± = — etc. 
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APPENDIX C 



We outline a rough estimate of the determinant of the fluctuation opera- 
tor (3.3), regulated by an IR cutoff L and a UV mode cutoff M. In effect, we 
restrict the eigenvalue range to [— |S +§] with open boundary conditions and con- 
vert to Fourier modes 

fa) = zY, e2nikx/L Pk (c.i) 

kez 

before setting p k = for \k\ > M. 

The Fourier transform of S'^(x,y) reads 

V -I0k+l,0 N [Pk+l - PkPl\ J 

which, for the simple case of V = \x 2 and after reinstating f(z) — > \nz 2 , can be 
evaluated explicitly, 



(C.3) 



- L t (4nk\ 

Pk = 2^fc J H — j 

fkl = ^(^ k ) + f) S k+l,0 = fkh+1,0 

with si(x) = — f x By some algebraic identities we rewrite 

det'S" = det' [5k+i,o-Nf k (p k+ i-p k pi)] = det[(5 k+ i fi -5 kfi 5i fi )-Ny^p k+ i^/Ji] 

(C.4) 

where the prime means omission of the zero mode (k — 0), and we have used po — 1. 
For large values of \k\, one has 

Pk - \f&W 2 ^{^-^) , f k « -jij- (C.5) 

so that the matrix elements for (C.4) fall off as |AH-/|~ 5 / 2 perpendicular to the 
diagonal. We can make a rough estimate of det ' S" by boldly taking into account 
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only the diagonal elements, i. e. 

« -iV/oIIa + f) 2 = 2iV( M ^) 2 (C.6) 

reg 



M 

det'5"' 

fc=l 



since fo = —2. Let us now choose a UV cutoff proportional to iV, say M = aN+f3. 
Via Stirling's formula then 



itr'lnS" 



« [(a+l)ln(a+l) -alna] iV + £>(lniV) . (C.7) 

reg 



More accurate estimates have to consider some sub-diagonals as well. 



2G 



REFERENCES 

1. M.L. Mehta: Random matrices, Academic Press (New York and London 
1967). 

2. E. Brezin and V.A. Kazakov, Phys. Lett. B236 (1990), 144; 
M.R. Douglas and S.H. Shenker, Nucl. Phys. B335 (1990), 635; 

D. J. Gross and A. A. Migdal, Phys. Rev. Lett. 64 (1990), 127,717; 

E. Brezin, M.R. Douglas, V.A. Kazakov and S.H. Shenker, Phys. Lett. B237 
(1990), 43; 

C. Crnkovic, P. Ginsparg and G. Moore, Phys. Lett. B237 (1990), 196; 
M.R. Douglas, Phys. Lett. B238 (1990), 176. 

3. E. Brezin, C. Itzykson, G. Parisi and J.B. Zuber, Comm. Math. Phys. 59 
(1978), 35. 

4. D. Bessis, C. Itzykson and J.B. Zuber, Adv. Appl. Math. 1 (1980), 109. 

5. B. Sakita, Quantum Theory of Many-Variable Systems and Fields , World 
Scientific (Singapore 1985). 

6. A. Jevicki, Nucl. Phys. B146 (1978), 77. 

7. A. Jevicki and B. Sakita, Nucl. Phys. B185 (1981), 89. 

8. J.D. Colin and S.P DeAlwis, IAS preprint IASSNS-HEP-91/7 (Feb. 1991). 

9. D. Karabali and S. Sakita, Int. J. Mod. Phys A6 (1991), 5079. 

10. S. Ben-Menahem, SLAC preprint SLAC-PUB-5377 (Nov. 1990). 

11. H. Sohngen, Math. Z. 45 (1939), 245. 

12. F.G. Tricomi, Quart. J. Math. Oxford Ser. (2) 2 (1951), 199. 

13. E.C. Titchmarsh, Introduction to the theory of Fourier integrals (Oxford, 
2nd ed. 1948) Chap. V. 

14. O. Lechtenfeld, IAS preprint IASSNS-HEP-91/2 (Jan. 1990), to appear in 
Journ. Mod. Phys. A. 

27 



